Self-consistent size and velocity distributions 
of coUisional cascades 



Margaret Pan^ and Hilke E. Schlichting^''^''^ 
mpan@astro . berkeley . edu , hilkeSucla . edu 

ABSTRACT 

The standard theoretical treatment of colhsional cascades derives a steady- 
state size distribution assuming a single constant velocity dispersion for all bodies 
regardless of size. Here we relax this assumption and solve self- consistently for the 
bodies' steady-state size and size- dependent velocity distributions. Specifically, 
we account for viscous stirring, dynamical friction, and collisional damping of 
the bodies' random velocities in addition to the mass conservation requirement 
typically applied to find the size distribution in a steady-state cascade. The 
resulting size distributions are significantly steeper than those derived without 
velocity evolution. For example, accounting self-con sistently fo r the v elocities can 



change the standard q = 3.5 power-law index of the iDohnanyil ( 1l969l ) differential 
size spectrum to an index as large as g = 4. Similarly, for bodies held together 
by their own gravity , the corresponding power-law index range 2.88 < g < 3.14 



of iPan &: Saril (120051 ) can steepen to values as large as g = 3.26. Our velocity 
results allow quantitative predictions of the bodies' scale heights as a function of 
size. Together with our predictions, observations of the scale heights for different 
sized bodies for the Kuiper belt, the asteroid belt, and extrasolar debris disks 
may constrain the mass and number of large bodies stirring the cascade as well 
as the colliding bodies' internal strengths. 



1. Introduction 

Believed to be a primary mechanism operating in circumstellar dusty debris disks as 
well as our own Kuiper and asteroid belts, collisional cascades — the transfer of mass from 
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larger to smaller sized bodies via collisions between those bodies — are ubiquitous in our 
galaxy. Their widespread occurrence and potential importance in understanding planet for- 
mation and planet-disk interactions have naturally provoked considerable study. Theoretical 
treatments predicting the collisional size distributioi i have incl u ded a nalytic work as well as 
numerical simulations. The pioneering treatment of iDohnanyil ( 1l969l ). who analytically cal- 
culated the size distribution for a steady-state cascade of constant-strength bodies, has been 
elaborated upon, extended to size-depen dent strength laws, and app li ed to different phys i- 



cal contexts by several authors, including IWillianis fc Wetherilll (ll994f):lTanaka et al.l % 
O'Brien fc GreenbereJ j2003l ): lKenvon fc BromlevlJ2004l ) : [OTrien fc GreenbereJ j2005l ): lPan fc Sari 
( 2005 ): Lohne et al. ( 20081 ) . Some of these also considered non-p ower-law features in the size 



distribution such as wa ves due to the gravity-strength transition (lO'Brien &: Greenberg 



2003 



20051 : iPan &: Sarill2005l ) or changes in the fragment size spectrum (IBelyaev &: Rafikovll201ll ). 
Many numerical studies of collisional cascades have also been performed (see , for example, 
Campo Bagatin et al. 1993; b urda fc Dermott 1997 Kenyon &: Bromlevll2004i iKrivov et al.l 



2005 



Kenyon fc Bromley ,2008i : iLohne et al 



2008 



Wvattll2008l : lFraser k Kavelaard l2009[ ) . 



The analytic and most of the numerical works on collisional cascades have generally as- 
sumed that the bodies' velocity dispersion is independent of size and constant in time once 
the cascade begins. Nevertheless, we expect processes like viscous stirring and collisional 
damping to affect all bodies' velocities throughout the cascade's lifetime. Here we incorpo- 
rate velocity evolution processes into our treatment of collisional cascades and find the size 
spectrum and size-dependent velocity dispersion self-consistently. In §2] we give expressions 
for the physical processes operating in the cascade. The well-known m ass conservation con- 
dition is the basis of previous work beginning with IDohnanyil ( 1l969l ). so we summarize it 
quickly in §2.11 In §2.21 we describe the stirring and damping processes affecting the veloc- 
ities, and in §2.31 we give expressions for the velocity equilibrium required in a steady state 
cascade. These are the velocity analogues of the mass conservation conditions of §2.11 In §3] 
we impose mass conservation and velocity equilibrium together to find the size and velocity 
power-laws of steady-state cascades. As we explain, in a disk of bodies with a power-law 
size distribution we expect to see up to three different velocity regimes. We derive velocity 
and size power laws in all three regimes for both gravity-dominated and strength-dominated 
bodies. In §5] we compare our analytic results to those of our fragmentation simulations. 
Finally, in ^ we summarize and discuss our findings. 
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2. Size and velocity evolution processes 

In order to find the size distribution and velocity function self-consistently, we assume 
a debris disk that occupies an annulus with typical orbital angular frequency fl. The bodies 
in the cascade have uniform composition and body mass density p, and their sizes r cover 
the range [rmin, ''^max]- We write the differential body size spectrum as dN/dr oc r~'^ and the 
velocity function as v{r) oc r^. We consider in turn how mass conservation in the cascade 
and velocity evolution via gravitational stirring and coUisional damping constrain q and p. 
Since our primary goal is to clearly delineate the relevant physical processes, we work to 
order of magnitude throughout. 



2.1. Mass conservation 

We begin with mass conservation, the basis for most analytic cascade treatments to 



date. Our discussion of mass conservation parallels that of iPan fc Saril (120051 ). and we refer 
readers to that work for a more detailed description. In a steady-state collisional cascade 
where mass is conserved in catastrophic collisions, the total mass destroyed per unit time 
per logarithmic interval in radius must be independent of size. This implies 

constant = pr^ ■ N{r) ■ \^ " ■ ■ v^^x (1) 

volume 

= pr" ■ N{r) ■ ^ ^ " (2) 

where rsir) is the size of the smallest body, or bullet, that can destroy a target of size r 
in a collision and v^^i is the typical relative velocity of bullets and targets. The second line 
follows because we assume isotropic velocities, so that the scale height of the disk is of order 
frei/n. The volume of Eq. [1] in which the bodies move is the area occupied by the disk 
midplane times this scale height, so the mass conservation relation for a disk is independent 
of velocity. 

We further assume that the way the bodies break is independent of size — that is, that 
the shape of the average fragment size distribution is size-independent. We parameterize the 
bullet-target size relation as a power law rsir) oc r". Then Eq. [2] yields 

6 + a 

? = -r— ■ 3 

1 + a 

The value of a depends on how much energy is lost in the post-impact destruction process. 
In the gravity regime, we can think of the destruction as a shock induced in the target by the 
bullet which propagates to the antipode of the impact site. The limiting cases for the shock 
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propagation are energy conservation and momentum conservation in the shocked material; 
these give respectively 



.3^2 
v,3 



pr^Vesc{r) ~ pr%v{r) 



a = (5 - 2p)/3 
a = {A-p)/3 



(4) 
(5) 



Here we assumed p > 0: p < would in principle require arbitrarily large velocities for 
arbitrarily small sizes, so we will not consider that case here. Numerical simulations of 
catastr ophic collisions find 1.37 < a < 1.66 with constan t collision velocities (see, for ex- 
ample, IStewart fc Leinhardtl l2009t iBenz fc Asphaud Il999l . and references therein); this is 
consistent with the range 4/3 < a < 5/3 for p = which we find above. Together, Eqs. |3HS] 
imply 

22 -p 23 - 2p 



> q > 
7-p 8-2p 

The inequalities hold if p < 1, which as we will see in % 



(6) 



lis satisfied. 



In the strength regime, a depends on the material properties of the body, which are 
often parameterized as Q*{r), the energy per unit mass needed to destroy a body of size r. 
In the strength regime sinaulatio ns find Q*{r) oc where 0>7>— 1/2 (IBenz &: Asphaug 
I999I : Istewart fc Leinhardtl l2009h . Then 



pr^Q*{r) = pr^v'^{r) — 

where, again, we assume p > 0. With Eq. [3] this gives 

21 + 7 - 2p 
6 + 7 — 2p 



a 



1 + 



7 - 2p 



(7) 



q 



As an example, iDohnanyil ( 119691 ) used in effect 7 
immediately yield g = 7/2 in Eq. [HI 



p 



in the strength regime; these 



2.2. Velocity evolution processes 

We now consider velocity evolution. Physically, v{r) depends on stirring from larger 
bodies and damping from collisions with and dynamical friction from smaller bodies. In the 
following, we explore the stirring-damping balance in detail. Motivated by observations of 
the asteroid and Kuiper belts, we work in the regime where the typical relative velocity in 
an encounter between two bodies is larger than either body's Hill velocity, which is of order 
a body's escape velocity times the one-sixth power of the ratio of the body's mass and the 
central mass. 
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We begin by writing expressions for the rates at which viscous stirring, colhsional damp- 
ing, and dynamical friction damping affect a body of size r. The rate at which bodies are 
viscously stirred by bodies of size R > r depends on which of them is moving faster. If 

Vesc{R) > v{R) > v{r), the rate is 



1 dv{r) 



v{r) dt 



N{R)R^n fv,,,{R)y fv,,,{R) 



area \ "^(-R) / V '^('") 



(9) 



This focusing factor applies because we need only double f(r), not necessarily v{R). If 
v{R) < v{r) < Vesc{R), we have 



1 dv{r) 



v{r) dt 



N{R)R^n fv,,c{R) 
Stir area V v{r) 



(10) 



Because Vesc{R) oc R, the largest bodies in the system do most of the stirring unless g > 5 
when v{R) > v{r) or unless q > 7 when v{R) < v{r). We expect these conditions to hold in 
real systems, so we will assume them throughout. We will show in ^that this assumption 
is self-consistent in the cascade. 

The rate at which bodies of size r are damped by direct collisions with bodies of size 
s < r is 



1 dv{r) 



v{r) dt 



damp area 



Nis)ns' . (11) 



There is no focusing factor here because, as we discuss later, dynamical friction damping 
is faster than damping by direct collisions only if v{r) < fesc('") 0- Note that bodies in 
the cascade must have v{r) > Vesc{r) because catastrophic collisions would otherwise be 
impossible. Eq. [11] implies that if g > 4, colhsional damping is dominated by the smallest 
bodies in the disk, giving s = rmin- 

When g < 4, collisions between equal-sized bodies dominat^, and Eq. [TT] becomes 



1 dv{r) 



v{r) dt 



N{r)Q 2 



damp area 



. (12) 



However, in the cascade rB{r) < r. Then the colhsional destruction rate will be at least as 
fast as the damping rate of Eq. [T2]as long as g > 1; the rates are equal only when rsir) = r. 



^Indeed, if v{r) > Vcscii^), dynamical friction is equivalent to elastic direct collisions. 

^The g = 4 size spectrum is a marginal case in which collisions with bodies of all sizes should contribute 
equally to the damping. Since this represents only an order unity correction to the damping rate, the scalings 
given remain valid. 
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Simulations of catastrophic collision ejecta indicate that nearly all the kinetic energy relative 
to the bullet-targ et center of mass is lost to heat in a catastrophic collision (see, for example, 
Jutzi et al.ll2010l ). Then bodies whose bullet-target mass ratio is not too small should lose 
most of their velocity dispersion in a catastrophic collision. If we assume that bodies are 
indeed damped whenever they are destroyed, then 



1 dv{r) 



v{r) dt 



1 dv{r) 



damp 



v(r) dt 



N{rB{r)yn 



coll 



area 



(13) 



If instead destructive collisions cannot damp effectively — that is, if the largest fragment 
retains most of its pre-collision velocity — then whether or not coUisional cooling is effective 
depends on the age of the disk. Bodies of size r are collisionally cooled only if the disk age 
is longer than the collision timescale implied by Eq. 



The dynamical friction damping rate of size r bodies by size s < r bodies also depends 
on whether v{r) > v{s) or vice versa. By analogy to the expressions of Eqs. P and [TUl for 
viscous stirring, the two expressions for dynamical friction are 



' N{s)ns^ 



1 dvir) 



v{r) dt 



area r 



df 



^^csc('^)y 


f Vcsc{r) 


v{s) ) 


\ v{r) 


N{s)VLs^ 


( Vesc{r) 


area r 


[ v(s) 



v{s) < v{r) < Vesc{r) 



(14) 



As mentioned above, dynamical friction acts faster than direct collisions by fesc(^) / (f (s)f (r))^ 
if v{r) > v{s) or by (vesd^) /v^s))* if f (r) < v{s), so it applies to bodies with t>(r) < fesc('^) 
which have not entered the cascade. The dynamical friction damping rate scales as s'^~''~2p jf 
v{r) > v{s) and as 54-9-4? ^^^.^ ^ v{s). Then the smallest bodies with s = rmin dominate 
the damping if g + 2p > 4 when v{r) > v{s) or if g + 4p > 4 when v{r) < v{s). 



2.3. Velocity equilibrium 



With expressions in hand for rates of velocity evolution, we can impose the steady-state 
condition that stirring and damping balance. In addition to the size and velocity power laws 
q, p of bodies in the cascade, we consider the analogous power laws q', p' for any bodies of 
size r > rmax which may be present in the disk but are too large to have entered the cascade. 
If the bodies in the disk formed through core accretion, we would ex pect 1 < g' < 5 (see, 
for example, iKenyon fc Bromley! l2004l . l2008l : ISchlichting fc Sarill201ll ) as well as the g < 5 
we already assumed. Motivated by the discussion after Eq. [101 we let the largest bodies in 
the disk have size R. 
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We first consider bodies in the cascade. As explained in §2.2[ these bodies are vis- 
cously stirred and collisionally damped. For cascades in which catastrophic collisions damp 
velocities effectively, velocity equilibrium means 







1 dv{r) 
v{r) dt 
1 dv{r) 



1 dv{r) 



v{r) dt 
' N{R)R'^Q 



v{r) dt 



(15) 



coll 



area 



v{r) 
N{R)R'^n 



area 



viR) 

VescjR) 

v{r) 



area 



V{r) < V{R) < Vesc{R) 



N{rB{r)yn 



area 



(i?) > v{r) > v{R) 



(16) 



where we have applied Eqs. IHl [IHl and[T31 Equivalently, if v{R) > v{r) the ratio of stirring 
and collision rates is 



N{R) 



N{rB{r)) \r 



R 



a(l-g)+2+2p 



(i?)V fVesciR) 
i?2 



v{R) 
NiR) 



A^(rB(rmax))rLx V HR) 



(17) 
(18) 



and if v{R) < v{r) this ratio is 

,a(l-g)+2+4p 



N{R) 



R^ 



VescjR) 
'"max \^(^max) 



(19) 

Note that we have transferred the coefficient of r"^^~'^^ in N{rB{r)) to A^(rB(rmax)) in the last 
step. In Eqs. [18] and [191 the second through last terms on the right-hand side are simply the 
ratio of the stirring and collision rates for size Tmax bodies, those at the top of the cascade: 
compare them, for example, to the right-hand side of Eq. [TTJ which is the ratio of stirring 
and collision rates for size r bodies. This indicates that if the stirring and collision rates for 
size Tmax bodics balance — which we expect since these bodies have just entered the cascade 
— the rest of the cascade will also be in velocity equilibrium if 



1 + 



1 + 



2 + 2p 

a 

2 + Ap 
a 



v{R) > f (r) 
v{R) < v{r) 



(20) 
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In the gravity regime, using Eqs. H] and |5] together with Eq. |20] gives 

10 + 5p ll + 4p 

^-->,>^-^ viR)>vir) 

10 + lip 11 + iop , • 

— >q> v{R) < v{r) 

A — p — 2p 

The inequahties for q hold when —1/2 < p < 1, which as we will see in ^ is satisfied. 
Similarly, in the strength regime we use Eq. [7] with Eq. |2D] to get 

9 + 7 + 4p 



3 + 7 - 2p 



v{R) > v{r) 

(22) 



q = — — v{R) < v{r) 

3 + 7 - 2p 

We next consider a disk in which catastrophic collisions do not damp the velocities. 
This may occur, for example, if rB{r) ^ r, in which case the center of mass velocity of 
a colliding bullet-target pair is dominated by the target velocity. Then conservation of 
momentum dictates that the velocity of the largest coUisional fragment will be quite similar 
to the target's velocity even if all of the relative kinetic energy between the bullet and target 
is lost. If g < 4 and if the system's lifetime is at least as long the timescale for two bodies of 
size Tinax to Collide, damping occurs through collisions between equal-sized bodies according 
to Eq. [121 This damping mechanism dominates for all bodies with v{r) > fesc(^)- While this 
condition holds over the entire cascade, it may hold for bodies outside the cascade as well. To 
see that all bodies in the cascade are included, note that if f (r) < Vesc{r), then the impact 
energy in a collision between equal-sized bodies, r^pr^v'^{r), is less than the gravitational 
binding energy ~pr^Uggp(r) of either body. If p > 0, the impact energy in a collision with a 
smaller bullet is likewise less than the gravitational binding energy of the target. Since both 
gravity-dominated and strength-dominated bodies require impact energy at least as large 
as their gravitational binding energies, v{r) > Vesc{r) is required in the cascade. However, 
f (r) > Vesc{r) may also apply for some bodies larger than rmax- 

A calculation entirely analogous to that of Eqs. [T61 - I20] above which uses damping by 
Eq. [12] rather than Eq. [121 and the size where v{r) ~ Vcsd^) instead than rmax, gives 

g = 3 + 2» v(R) > v(r) 

\ J \ J ^ (23) 

q = 3 + Ap '^(-R) < ^{r) 

If g > 4, we substitute Eq.[TTlfor Eq. [T^lin the above calculation. In this case N{r) disappears 
from the ratio of stirring and damping rates and we get a condition on p alone: 

p=l/2 v(R)>v(r) 

^ ' ^ ' ^ . (24) 

p = 1/4 v{R) < v{r) 
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If the cascade lifetime is short compared to the timescale for colhsions between bodies of 
size Tmax; then some bodies near the top of the cascade will not have had time to damp. For 
these undamped bodies, we expect shallower velocity power laws. 

Finally, we consider any bodies in the disk whose velocities are smaller than their escape 
velocities. They cannot be part of the cascade, so we expect them to have sizes r > rmax- 
The equilibrium velocities of these bodies follow from a balance between viscous stirring and 
dynamical friction. If v{r) > v{R), we equate the stirring rate of Eq. [Hlwith the damping 
rates of Eq. [H] to get 

v{s) < v{r) 

(25) 

v{s) > v{r) 

Here s is the size of bodies which dominate the dynamical friction. Because we have broken 
power-law size and velocity distributions, and because the power-law breaks do not occur at 
r, we expect s to be independent of r. If v(r) oc r^', Eq. [251 implies 

p' = -3/2 v{s) < v{r) 
p' = —3/4 v{s) > v{r) 

This is indeed consistent with the v{r) > v{R) we assumed. Note that the kinetic energy 
per body, r^pr^v'^{r) oc r^"*"^^', cannot increase with decreasing body size, so these bodies 
lie outside the cascade. For the same reason we can neglect any dynamical friction heating 
effects, which contribute at most an order unity correction. 

If instead we assume v{r) < v{R) and replace Eq. [9]with Eq. [TO]above, no self-consistent 
solution for p' is possible. 

3. Steady-state size and velocity distributions 

We now solve simultaneously the mass conservation and velocity stirring/damping bal- 
ance conditions of ^ to find the steady-state size and velocity distributions in the disk. 
We first confirm that the steady-state condition — equivalent to requiring that A^(rmax) 
changes on a timescale long compared to collisions between and stirring of smaller bodies 
— is physical. Since the stirring cross-section of size r bodies scales as r~^^, smaller bodies 
are indeed stirred faster than the largest bodies break. Similarly, since smaller bodies have 
more total surface area than larger bodies as long as g > 3, smaller bodies break faster than 



v{r) / N{R) y^ 

v{r) f N{R) V'^ R^/^ 
^ ^ ( N{s) ) s3/4r3/4 



(26) 
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larger ones. Our assumption of a steady state is therefore reasonable for all bodies smaller 
than Tinax- Said another way, r^ax corresponds to the location of t he break seen in coUi- 



sional size distributions separa t ing; coUisional ari d primordial bodies (jO'Brien fc Greenberg 



2OO3I : iKenvon fc BromlevlEooi : |Pan fc Saril l2005l ). As bodies of size rmax break and N{r 



decreases, the normalization of the cascade below rmax should follow adiabatically. 

In addition, we can see that these solutions are stable by considering a perturbed cas- 
cade. If for any reason stirring becomes faster than catastrophic collisions, the velocities will 
increase and rs(r) will decrease until the catastrophic collision rate equals the new stirring 
rate. Similarly, if stirring becomes slower than catastrophic collisions, the velocities will slow 
and rB(r) will increase until collisions just balance stirring as long as v{r) > Vesd'^')- The 
timescale for bodies smaller than some size r to relax to this solution should be of order a 
few catastrophic collision times for size r bodies. 

We frame our discussion of the solutions via the velocity stirring/damping equilibria 
listed in §2.31 They suggest that given a disk in which a single rsir) power-law relation 
applies to all bodies, and in which cooling has had time to operate, up to three velocity 
regimes occui§|. First, the largest bodies, which are stirred viscously and damped by dy- 
namical friction, have velocities that are below their escape velocities but that increase with 
decreasing size according to Eq. [261 At the size for which the bodies' velocity equals their es- 
cape velocity, dynamical friction can no longer cool efficiently and the second regime begins. 
Bodies slightly smaller than this first transition size have velocities faster than both their 
own escape velocities and the velocity of the largest bodies in the disk. In this regime stir- 
ring proceeds according to Eq. [10] and damping proceeds by collisions. Third, if the cascade 
includes sufficiently small bodies, we expect for p > that the smallest bodies' velocities fall 
below the velocity of the largest body. In this regime stirring proceeds according to Eq. [9] 
and coUisional damping continues. While we would expect the power-law breaks associated 
with transitions between regimes will produce waves in the size and velocity distributions, 
we also expect that, on average, the sizes and velocities in each regime will be consistent with 
the q and p values we find. We discuss these three regimes — first, v{r) < Vesc{r)] second, 
f (r) > Vcsc{r) and v{R) < f (r); third, f (r) > Vesd^) and v{R) > v{r) — in turn below. 

The regime containing the largest bodies of sizes r > r^ax is simplest. Because its bodies 
are not part of the cascade, we cannot constrain their size distribution except to require that 
it satisfy the conditions for Eq. [10] to hold. Instead we expect that their size distribution 



•^If significant external stirring has occurred, not all of these three regimes may occur. For example, if all 
the bodies in the disk have velocities larger than their own escape velocities, dynamical friction will never 
be important. 
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N{r) oc has not changed since their formation. Regardless of whether the bodies are 
gravity- or strength-dominated, Eq. [26] gives 

p = -3/2 , l<q<7 v{s) < v{r) ^^^^ 

p = -3/4 , 1 < g < 7 v{s) > v{r) 

The remaining two velocity regimes, v{r) < v{R) and v{r) > v{R), may support cas- 
cades. We consider cascades with four different categories of rsir) relations characterized by 
1) whether the bodies are gravity- or strength-dominated and 2) whether the bullet-target 
size ratio is close enough to unity for catastrophic collisions to provide effective cooling. 

First we assume cooling by catastrophic collisions. This case requires a cascade. In the 
gravity regime, we solve Eqs. M and [20] simultaneously to get 



17-^241 ^241- 9 , 4-p 

p = , q = for a = 

^ 4 ' ^ 2 3 

11-^85 785-3 , 5-2p 

p = , q = for a = — - — 



31 7865 - 23 , 4-p 

p = , q = for a = 

^8^2 3 

1 . r 5-2p 

p = - , q = o for a = — - — 



v{R) > v{r) {2h 



v{R) < v{r) . (29) 



4 



This implies 



0.37 <p< 0.45 
3.26 > q> 3.11 

0.20 < p < 1/4 
3.21 > g > 3 



v{R) > v{r) (30) 



v{R) < v{r) . (31) 



In the strength regime, we likewise solve Eqs. [8] and [20] together for 

9 + 7-^69 + 67 + 72 -1 - 7 + J69 + 67 + 72 
p = , q = v{R) > v{r) 



(32) 



15 + 27- v/201 + 367 + 472 -7 - 27 + v/201 + 367 + 472 
p = , q = 'v{R) < v[r) 



(33) 
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For the range — 1/2 < 7 < 0, these give 

0.090 <p< 0.17 . ^ . ^ 

v(R) > v(r) 34 

3.82 >q> 3.65 
0.054 <p< 0.10 

v(R) < v(r) . (35) 

3.78 >q> 3.59 

Now we assume catastrophic colhsions cannot damp the velocities significantly, so that 
the cooling timescale is the time it takes for a given body to collide with a total mass equal 
to its own. We also assume the lifetime of the disk is longer than this cooling timescale for 
all bodies with v{r) > Vcsc- In the gravity regime, these bodies all participate in the cascade: 
''^(^max) ~ ^cscl'^max)- Their Steady state sizes and velocities should follow from Eqs. [H] and 
When a = {A — p)/3 this gives 



(36) 



p = , q = 9 — y/SA for a = — — "^iR) > v{r) 



'34 4 — p 
p = , g = 9 — Vm for a = "^(-R) < '^{f) 



We find, however, that no solution with p > and g < 4 is possible when a = (5 — 2p)/3. It 
turns out g > 4 is also impossible in the gravity regime. If g > 4, Eq. [6] implies p > 3/2, and 
since fesc('^) oc r, having p > 1 means that v{r) will fall below Vcsdr) at some r, stopping 
the cascade. Then the maximum a allowed must lie between (4 — p)/3 and (5 — 2p)/3. To 
find this limiting value, we recast Eqs. |l]and|5]as 

pr'^vLir) - prlv^ir) ^ a = 1 + ^^(l - p)/3 (37) 
where 1 < (3 < 2. With Eqs. [3]and|23l this gives 



3 ^ v/36-6/3 + 4/32 

^ " ^ 2P '"^^ '"^^^ 

12 + /3 ^ v/l44 + 12/3 + 9/32 

P = ± v{R) < v{r) 



(38) 



A look at the zeros of dp/dP shows that p is monotonic for the relevant /3, so the limiting 
a and (3 should occur at a limiting value of p. For gravity-dominated bodies, < j9 < 1 as 
discussed above; for both v{R) > v{r) and v{R) < v{r), the only /3 between 1 and 2 that 
satisfies p = 0orp=lis/3 = 3/2 at p = 0. When p = 0, Eq. [231 gives g = 3. Then the 
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allowed p, q in the gravity regime are 



< p < 0.085 
3 < g < 3.17 
< p < 0.042 
3 < g < 3.17 



v{R) > v{r) 



v{R) < v{r) 



(39) 



(40) 



In the strength regime, not all of the bodies with v{r) > Vesc{f) can participate in 
the cascade. For those in the cascade, we first assume g < 4 and solve Eqs. |8] and 
simultaneously. This gives 



P 



P 



4 + 7-^/4+167 + 72 



10 + 7 - ^4 + I67 + 



r 



v{R) > v{r) (41) 



^±2 \/l9 + 147 + 7 _2 ^ ^ ^ ^ ^ ^ ^ (42) 



When 7 < 0, the only real solutions to Eq. HUhave g > 4, which is inconsistent. The allowed 
ranges in p, q are 



p = l/2 
g = 4 



v{R) > v(r) 
7 = 



(43) 

(44) 



1/4 >p> 0.16 v{R)<v{r) (45) 

4 > g > 3.64 -1/2 < 7 < . (46) 

The g > 4 which arose above when v{R) > v{r) and 7 < suggests that we look for a 
solution where the smallest bodies in the system dominate the collisional damping — that 
is, a solution using Eq. |2l] instead of Eq. [231 Indeed, Eqs. [8] and |24] together give 

20 + 7 

p = 1/2 , g = — -2 v{R) > v{r) (47) 



5 + 7 



and, for -1/2 < 7 < 0, 



p=l/2 
13/3 > g > 4 



v{R) > v{r) . (48) 



For bodies with v{r) > Vesd^) but r > rmax — those not in the cascade — the primordial 
size distribution g' applies. As long as g' satisfies the conditions on Eq. [9l the velocities follow 
from this and Eq. [23] if g' < 4 or Eq. [2l if g' > 4. 
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Finally, if the coUisional cooling timescale is shorter than the age of the cascade, the 
velocity distribution will be shallower than predicted in the relevant regime above. How much 
shallower depends on particulars of the stirring timescale and the energy loss per collision. 
For example, if the kinetic energy lost in a catastrophic collision is so small that the kinetic 
energy retained by the fragments is larger than the energy they gain via stirring in one 
collision time, p will instead depend on exactly how much energy is lost in an average collision. 
In turn the energy loss per collision depends hea vily on the bodies' internal structure, which 
is very poorly constrained (ILeinhardt et al.ll2008l . and references therein). We will not discuss 
this uncooled regime in detail here. 

Our results for p and q in all the velocity and strength law regimes discussed in this 
work are summarized in Table [T] Note that all the size distributions are steeper than those 
that obt ain when fixed ve locities are used (p = 0); these are 3.14 > g > 2.88 for the gravity 
regime ( jPan fc Saril 120051 ) and 3.72 > q > 3.5 for the strength regime. The steepening is 
certainly consistent with smaller velocities for smaller bodies: lower velocities mean larger 
bullets are needed to break a target of a given size; an increase in bullet size corresponds 
to a decrease in the number of bullets for q > 0; and a steeper size distribution offsets this 
decrease. Tabled] also confirms that our assumption p < 1 of §2. H is self-consistent. 



4. Comparison with numerical simulations 

To test the analyti c results above we used a numerical cascade simulation based on 



the coagulation code of ISchlichting fc Saril (120111 ). Because our goal here is to study the 



dominant physical processes in the cascade — viscous stirring, coUisional and dynamical 
friction damping, and mass transfer from larger to smaller body sizes — we neglect factors 
of order unity in the stirring and damping rates. We study a single belt of bodies orbiting 
in an annulus about a much more massive star. We take the initial total mass in bodies to 
be about IMsarth, and we assume the bodies have bulk density 1 g/cc and follow the mass 
and velocity evolution of bodies with radii ranging from 1 m to 3000 km, a few times the 
size of Pluto. 

As a first test of our velocity evolution theory, we artificially fix the size spectrum in the 
simulations and allow only the velocities to evolve. In Figure [1] we show as an example the 
test results with strength-dominated 7 = bodies and a fixed q = 3.6 size spectrum. Since 
we do not allow for catastrophic collisions in this run and since we fix the size spectrum at 
g < 4, the coUisional damping is dominated by collisions between similarly sized bodies as 
given in Eq. [121 The resulting steady-state velocities obey Eq. [231 which for q = 3.6 means 
p = 0.3 if v{R) > f (r) and p = 0.15 if v{R) < v{r). We expect the velocities for large bodies 
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damping 
mechanism 


v(R) > v[r) 


v[R) < v{r) 


reierences 


v{r) > Vesc{r): 
includes all 


gravity 
regime 


catastrophic 
collisions 


0.37 <p< 0.45 
3.26 > q> 3.11 


0.20 <p < 1/4 
3.21 > q>3 


Eqs. M M 


bodies in 
cascade 




collisions with 
equal-sized 
bodies 


< p < 0.085 
3<q< 3.17 


0<p< 0.042 
3<q< 3.17 


Eqs. [391 M 




strength 
regime 


catastrophic 
collisions 


0.090 <p< 0.17 
3.82 > q> 3.65 


0.054 <p< 0.10 
3.78 > q> 3.59 


Eqs. [Ml M 






collisions with 
equal-sized 
bodies 


p = l/2 
q = i 


1/4 > p > 0.16 
4>q> 3.64 


Eqs. [Ml M 






collisions with 
smallest bodies 


p=l/2 
13/3 > q>i 




Eq.M 


v{r) < Vescir): 
bodies too 
large for 
cascade 


gravity or 

strength 

regime 


dynamical 
friction^ 




p = -3/2 or p = -3/4 
1< g < 7 


Eq.m 



Table 1: Summary of velocity power laws p and size power laws q in steady-state for all of 
the stirring and damping regimes discussed in this work. 



°p — —3/2 applies when Vcsc{r) > v{r) > v{s); p = —3/4 applies when v{r) < v(s) < Vescir). Here v{s) is the 
velocity of the bodies providing the dynamical friction. 
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with velocities below their own escape velocities to follow Eq. |26l Our simulations agree well 
with these numbers. 

Similarly, we test our mass cascade implementation by artificially fixing the velocity as 
a function of size and allowing only the size distribution to evolve. Figure [2] shows the results 
of a test with strength-dominated 7 = bodies and velocities fixed to a broken power law 
with p=l/4, p=l/8. For these Eq. Ogives q = 3.73, q = 3.61; these agree well on average 
with our simulations. Our simulations also show waves as mentioned in §31 these are induced 
by the break in the velocity distribution as well as the artificial "breaks" in the mass power 
law created by the finite range of body sizes in our simulations. 

Finally, we allow both the size and velocity distributions to evolve in the simulations. 
Figure [3] shows an example again using strength-dominated bodies with 7 = 0. We assumed 
in this run that the coUisional damping of the velocity dispersion is dominated by collisions 
between like-sized bodies (see Eq. [12]). This criterion applies when catastrophic collisions 
do not damp the velocity dispersion significantly, which may occur for small bullet-to-target 
ratios. Here the steady-state solution of Eq. HU applies, and 7 = implies p = 1/2, q = 4 
when v{R) > v{r) and p = 0.16, q = 3.64 when v{R) < v{r). Again, these agree well with 
our simulations on average in each of the three different velocity regimes. 

This model and the results shown in Figure [3] may, for example, apply at the end of 
protoplanetary growth in a planetesimal disk. Initially, the velocity dispersion is so small that 
collisions lead to growth. As the largest bodies — "protoplanets" — grow, they continue 
to excite the small planetesimals' velocity dispersion; their velocities grow on the same 
timescale as the large protoplanets ' sizes (for a comprehensive description of this growth 



phase see ISchlichting fc Saril (120111 )). Once the system reaches an age comparable to the 
small planetesimals' collision time, but before collisions become destructive, the balance 
between gravitational stirring and collisional damping determines the planetesimals' velocity 
dispersion. This phase is similar to the situation shown in Figured! but with a mass spectrum 
that continues to evolve due to planetesimal accretion. Finally, the planetesimals' velocity 
dispersion is excited sufficiently above their escape velocities that destructive collisions set 
in. This stage is shown in Figure [31 The mass spectrum now no longer refiects the growth 
history; instead it is determined by the collisional evolution. 



5. Summary 



We have found self-consistent steady-state solutions for the velocity function and size 
distribution of collisional cascades in the super-Hill regime. These solutions occur when 
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Fig. 1. — Comparison between analytic results (dashed blue line) and simulations (solid 
orange curve) for steady-state velocities in a system of strength-dominated bodies with 7 = 0. 
The top panel shows the mass in a given log2 mass bin as a function of radius, which in this 
run is fixed with q = 3.6; the bottom panel shows the simulations and analytic results for the 
velocities. Since we do not allow for catastrophic collisions in this run and since we fix the 
size spectrum at g < 4, the coUisional damping is dominated by collisions between similarly 
sized bodies as given in Eq. [121 There is good agreement in each of three velocity regimes. 
The smallest bodies, which have velocities greater than their own escape velocities but less 
than v{R), follow p = 0.3 (see Eq. [23l) . Larger bodies still small enough to have velocities 
larger than their own escape velocities, but which have velocities greater than v{R), follow 
p = 0.15 (see Eq. [23]). Finally, the largest bodies have p' = —3/2 because they are subject 
to dynamical friction by small bodies with velocity dispersion v{r) < v{R) (see Eq. [26l) . 
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Fig. 2. — Comparison between analytic results (dashed blue line) and simulations (solid 
orange curve) in a steady-state system of strength-dominated bodies with 7 = and fixed 
velocity distribution. The top panel shows the resulting mass spectrum plotted as mass in 
a given log2 mass bin as a function of radius (solid orange curve) and the corresponding 
analytic results (dashed blue line); the bottom panel show the fixed velocity distribution. 
We begin the simulations with an initial size distribution q = 3.6, which evolved to the 
q = 3.61 for v{r) > v{R) and q = 3.73 for v{r) < v{R) as expected from Eq. [8]with p = 1/4, 
p = 1/8. Waves are clearly visible as oscillations in the steady-state mass spectrum. We note 
here that the wavelength of the waves changes as one transitions from the p = 1/4 to the 
p = 1/8 velocity distribution. This change in wavelength refiects the velocity dependence in 
the bullet-to-target ratio. 
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Fig. 3. — Comparison between analytic results (daslied blue line) and simulations (solid 
orange curve) in which the size distributions and velocities are both evolved together. The 
top panel shows the mass spectrum; the bottom panel shows the velocity distribution. In 
this run we assumed that the coUisional damping of the velocity dispersion is dominated 
by collisions between like-sized bodies (see Eq. [T2l) . This damping criterion applies when 
catastrophic collisions do not damp the velocity dispersion significantly, which may occur 
for small bullet-to-target ratios. The agreement between the simulations and our analytic 
results in Eq. HTJand Eq. |26]is good on average for the v{R) < v{r) and v{K) > v{r) regimes 
in both the mass and velocity plots as well as for the v{r) < Vcsc{i^) regime, (see caption 
of Figure [T] for a description of the regimes). Waves due to both the transition between 
v{R) > v{r) and v{R) < v{r) and the ends of our simulation range are again visible in the 
data. 
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mass conservation is satisfied and when viscous stirring balances velocity damping. Three 
kinds of velocity equilibrium may occur. For the biggest bodies, which have velocities slower 
than their escape velocities, viscous stirring and dynamical friction balance. These bodies' 
velocities increase with decreasing size until the size at which the velocity and escape velocity 
are equal. Since dynamical friction is inefficient for bodies with velocities faster than their 
escape velocities, stirring balances damping by direct collisions for all smaller bodies. Bodies 
just smaller than this first transition size have velocities faster than both their escape veloc- 
ities and the velocity of the largest bodies in the system. A second transition occurs at the 
body size whose velocity equals that of the largest bodies in the system. Bodies smaller than 
this second transition have velocities slower than the largest bodies in the system, so their 
stirring requires a different cross-section. The resulting size distributions for the gravity- 
and strength-dominated regimes are steeper than the ones expected with size-independent 
velocities. We find good agreement between the predictions of our theory and the results of 
our numerical simulations. 

To our knowledge, previous a nalytic t reatni ents of coUisional cascad e s hav e not consid- 
ered velocity stirring or damping. IWyattI ( 120081 ) and Kennedy fc WyattI ( I2OIOI ) study disks 
in which the cascade start time depends on orbital radius because the large bodies needed to 
excite the velocity dispersion and initiate a cascade take longer to accrete at larger orbital 
radii. However, they do not consider the effects of stirring or damping on colliding bodies' 
velocities as the cascade proceeds. 



Kenyon fc Bromleyl (120081 ) do account for simultaneous velocity and size spectrum evo- 



lution in their coagulation/fragmentation code. Our results here are not directly comparable 
to theirs because they do not account for energy lost during catastrophic disruptions and 
because the largest bodies in their simulations continue to accrete while their coUisional cas- 
cades operate. We plan to extend and modify our calculations to enable comparison with 
their findings. Other areas for future investigation include incorporating velocity stirring and 
damping into coUisional cascades covering both the gravity and strength regimes as well as 
the waves induced in the size and velocity power laws due to transitions between regimes. A 
good knowledge of the size and velocity distributions will also allow us to predict observables 
such as the dust production rate as a function of time and the scale height of the disk as a 
function of size or, for the smallest bodies, observing wavelength. 

Ongoing surveys of the Kuiper and asteroid belts provide observational size distribution 
and velocity data to which we can compare our results. Kuiper belt surveys indicate that 
its size spectr um follows a broken power law who se break falls at a body size of several tens 
of kilometers (IBernstein et al.l 120041 : iFraserl |2009| ): this break is interpreted as the top of a 
coUisional cascade. Typical Kuiper belt velocities are about 1 km/s, of order 30 times larger 
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than the escape velocities from the largest bodies in the cascade, so the typical bullet/target 
size ratio is far from unity. Then cooling by catastrophic collisions should be ineffective. 
Also, the timescales for the observed Kuiper belt objects (KBOs) to collide with bodies of 
equal size are longer than the age of the solar system. These KBOs have therefore not had 
time to cool; we would expect their average velocities should be very similar to those of the 
primordial KBOs. Indeed, small KBOs' eccentricities and inclinations show no significant 
trends with size. As for the size distribution, assuming the break exists, sur veys find a range 



of size distributions 1.9 < q < 3-9 for KBOs smaller than t he break size (IBernstein et al. 



2004J : iFraser et al.l l2008l : iFraserl l2009l : ISchhchting et al.ll2009f ). This is consistent with the 
2.88 < q < 3.14 we expect if p ^ but not strongly constraining. In the asteroid belt, 
typical relative velocities of ~5 km/s suggest catastrophic collisions are likewise ineffective 
at cooling. Surveys of the asteroid belt indicate a size distribution of q ^3.5 for large bodies 
of H magnitude smaller than about 15, or size larger than ~1 km (iGladman et al.l 12009 



and references therein). For smaller bodies, however, the slope become s shallower; different 
surveys report slopes ranging from q = 2 to q = 2.8 ( llvezic et al.l 1200 ll : lYoshida et al.ll2003l : 



Yoshida fc Nakamurall2007l : IWiegert et al.ll2007l ). While the overall size distribution slope is 
roughly consistent with the expected 2.88 < q < 3.14, we would predict that the average 
slope steepen for bodies smaller than about 100 m in size. Still, our theory alone suggests 
several possible causes for waves that might explain the observed break and its location. 
This again makes strong constraints difficult without further data on smaller bodies. 

We look forward to future observations of smaller KBOs and asteroids whose cooling 
time may be shorter than the belts' lifetimes and which will provide a longer size baseline 
with which to compare our theory. Future survey results of this kind will provide more 
stringent tests of our results and may shed light on the catastrophic collision process in our 
solar system. In particular, measurements of the slopes of the size and velocity distributions 
would provide a direct probe of the bodies' strengths. Similarly, observations of debris disk 
scale heights as a function of wavelength at millimeter wavelengths, for example with ALMA, 
would provide direct tests of our velocity power laws as well as constraints on the internal 
strengths of pebble-sized particles in those disks. 
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